Lifetime and fragment correlations for the two-neutron decay of °0 ground state 
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The structure and decay of ^®0 are studied in a three-body ^^0+n+n model. We have found 
extremely strong effect of the subbarrier configuration mixing on the width of true 2n emitters due 
to core recoil and neutron-neutron final state interaction. This effect is far exceeding analogous effect 
in the true 2p emitters. Our calculations provide reasonably narrow boundaries for the lifetime vs. 
decay energy dependence for the true 2n emission. An upper limit of ~ 1 keV for the decay energy 
of the unbound ^^O is inferred basing on the recent experimental lifetime value. 
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I. INTRODUCTION 

True few-body decay is an exclusive quantum mechan- 
ical phenomenon comprising simultaneous emission of 
several particles. Such a phenomenon leads to unusual 
lifetime systematics and complex correlation patterns for 
decay fragments. It takes place when the sequential emis- 
sion of particles is prohibited by specific energy condi- 
tions. These conditions appear to be widespread beyond 
the nuclear driplines. Two-proton (2p) radioactivity [l| 
and "democratic decay" are the best-known specific 
examples of a broader true few-body decay phenomenon 
which were under scrutinous investigation in the last 
decade Q. Nowadays, the studies of the neutron- rich 
systems beyond the proton dripline arc very active. Re- 
cently there is an outburst of experimental results con- 
cerning the light true two-neutron (2n) emitters: ^"^He 
M, ''Li ii, i6Be El, and ^^O ^MM- 

The possibility of the very long-living (radioactive) 
neutron emitters was considered in the framework of sim- 
plified model in our work [lj|. The true four- neutron 
emission was found to be the most prospective phe- 
nomenon for an experimental search. However, broad 
boundaries for existence of radioactive 2n emitters were 
also estimated, naming the ^^O system, information on 
which was absent at that moment, as one of the candi- 
dates. Since that time, upper limits were found experi- 
mentally for the decay energy of ^^O: Et — ISOl^g keV 
[ll| and Et < 120 keV 'l^. For known ^sq ground 
state decay energy of 770tio ^'^^ ^-'^\^ ^^^ ^^^ i^ clearly 
ascribed to be a true 2n emitter, see Fig. [T^a). 

Recently the halflife of ^^O was reported to be within 
the radioactivity timescale: T1/2 = 4.5l^;5(stat)±3(sys) 
ps [ij|. Taking into account experimental complexities 
and novelty of the method the Authors express cautious 
optimism about possibility of two-neutron radioactivity 
observation in their work. Prospects of the discovery of 
this new type of the radioactive decay call for further fo- 
cused experimental search and deeper theoretical insight, 
providing the guideline for such a search. This inspired 
us to perform more detailed studies of the long-living 
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FIG. 1. (a) The energy scheme for 2n decay of ^®0. (b) 
Schematic comparison of two decay mechanisms: (i) in two- 
body decay the angular momentum I is conserved in the pro- 
cess of decay; (ii) in three-body decay of the [d'^Ji^o initial 
configuration, the subbarrier configuration mixing boosts the 
penetration drastically. 



true 2n emitters by example of the ^^O system. Our 
work unveils complicated dynamical effects which should 
be accounted in calculations of true 2n emitter's widths. 
The important peculiarity of such a decay mechanism for 
^^O with [cP] structure is schematically illustrated in Fig. 
nib). 



II. THEORETICAL MODEL 

In studies of long-living true 2n emitters we use the 
three-body hyperspherical harmonics (HH) cluster model 
developed for studies of 2p radioactivity ^] . The model 
demonstrated applicability and high precision for broad 
range of de cay energies, masses, and structures of 2p pre- 
cursors 0, [la] ■ In this approach the following sequence 
of problems is solved: 

(iJ3-ST)*box = 0, (1) 

{H3 - Et)¥+^ = -i(rarb/2)*box , (2) 

T = j/N, (3) 

where the eigenenergy Et and corresponding eigenfunc- 
tion ^box are found as a solution of the homogenous 



Schrodinger equation ([T]) with some "box" boundary con- 
ditions; next, the wave function (WF) vl/(+) with a pure 
outgoing three-body asymptotic is found by solving the 
inhomogeneous Schrodinger equation ^ with an arbi- 
trary value of the width Farb! finally, the actual width 
F is defined in ([3]) via the outgoing fiux j and internal 
normalization N connected with WF ^(+). The width 
values obtained in model calculations are shown in Figure 

m 



III. THREE-BODY CORRELATIONS 

The three-body calculations in the HH method utilize 
the collective coordinates: the hyperradius p (describing 
collective radial motion) and the hyperangle 9p (respon- 
sible for geometry of the system at given p): 
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y Ai + A2 + .43 V Ai + A2 
p = v x^ + y'^ , dp — a.Tcta.n{x/y) . 

The boundary conditions of the HH method for the 
three-body decays without Coulomb interactions are well 
known. However, several effects, including the large scat- 
tering length in the n-n channel, lead to slow radial and 
basis convergence of the calculations. The hyperradii up 
to 2000 fm are used depending on the decay energy. Some 
examples of the basis convergence are given in Figure [31 
In the momentum space the three-body energy-angular 
correlations (which actually can be measured, see Q) are 
defined by the energy distribution parameter e and the 
angle 9k between Jacobi momenta k^: and kj, (conjugated 
to Jacobi radii x and y): 
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where M is an average nucleon mass. For systems con- 
sisting of two nucleons and "core" these correlations 
could be defined in two Jacobi systems: "T" (particle 
3 is the core) and "Y" (particle 3 is one of the nucleons). 
Coordinate and momentum space correlations for several 
cases of ^^O calculations are given in Figure |4l 



IV. POTENTIALS 

Taking into account the limited information on ^^O- 
n subsystem and exploratory character of our current 
studies, we use simplified potential sets. The n-n poten- 
tial, acting only in s-wave, is taken in the Gaussian form 
V{r) = -31exp[-(r/1.8)2]. For the "^^O-n channel we 



use a Woods-Saxon potential depending on the angular 
momentum {Is forces are not used). The only available 
data on ^^O spectrum [I4I allows to identify the ground 
state energy of 770+1° keV and the width of 172 ± 30 
keV which reliably indicates a d-wave state. Therefore, 
for the d-wave in the ^*0-n channel we use the poten- 
tial with depth Vd = —35.5 MeV, radius tq = 3.5 fm, 
and diffuseness a — 0.75 fm, which well reproduces the 
experimental properties of the resonance. The potential 
parameters for s- and p- waves used in different model cal- 
culations are given in Table HI The three-body potential 
depending on p with the typical Woods-Saxon parame- 
terization Q is used to control the decay energy Et in 
the systematic lifetime calculations. 



V. CALCULATIONS WITHOUT n-n FINAL 
STATE INTERACTION (FSI) 

The structure of ^^O is expected to be dominated by 
the [cP] configuration. If the decay proceeds also via the 
[(P] configuration then the three-body model should pro- 
vide the decay widths which have to be very close to 
the [(P] estimate results from |14|. To reproduce formal 
conditions of this estimate in the three-body model, two 
actions are required: (i) the n-n potential should be put 
to zero and (ii) the core mass should be infinite (actually 
the mass number ~ 100 is a sufficient approximation). 
One can see in Fig. [2] that such three-body calculations 
(the solid curve) are very close to the generic [(P] esti- 
mate, well following the energy trend F ~ Ej, expected 
for the [(P]q decay. 

The correlation patterns for this case are given in Fig. 
^a) . One can see that the "triple- ridge" correlation con- 
nected with [(P]q coupling in the internal region survives 
to the asymptotics producing very symmetric "triple- 
peak" correlation in the energy distribution between neu- 
trons in the Jacobi "T" system. 



VI. THE RECOIL EFFECT 

The surprise comes when the core recoil is treated cor- 
rectly with the right mass of the ^^O fragment used in 
the "no FSI" calculations. The accounting for the re- 
coil results in several orders of magnitude larger widths 
(the dashed curve in Fig. [5]). The slope of the curve also 
evolves with the decay energy Et ■ At the highest energy 
it is parallel to the [(P]o curve from ,14]. With energy 
decrease it changes first to [p^] then to s-wave slope at 
very small Et- This indicates that the core recoil leads 
to configuration mixing allowing the WF to "migrate" 
into partial wave configurations with lower centrifugal 
barriers thus boosting the penetration, see Fig. [Ijb). 

Partial wave analysis and also the analysis of correla- 
tions confirm this idea. One can can see in Fig. IDJb) that 
the triple-ridge [(P]o internal configuration is smearing 
out with radius increase and at about p = 150 fm it is 
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FIG. 2. (Color online) Width (halflife) of ^^O as a function of its decay energy in different model calculations. Gray curves show 
the estimates of Ref. ;141] for decay via pure orbital [i ] configurations coupled to the total angular momentum L = 0. Solid 
curve shows the "no FSI, no recoil" case, dashed curve corresponds to "no FSI" situation. Dotted curve shows the results with 
n-n FSI scaled by factor 0.25. Dash-dotted and dash-double-dotted curves correspond to "strong repulsion" and "moderate 
repulsion" in s- and p- waves in ^^O, while short-dash curve shows the case with "no repulsion in p-wave". Hatched areas give 
the experimental limits from Refs. [111 . \\3{ and the realistic theoretical limits from this work. 



replaced with a double-ridge configuration typical for a 
p-wave dominance. On asymptotics, this is reflected by 
a double-peak energy correlation between neutrons. 



VII. EFFECT OF THE n-n FSI 

It is found to be much stronger in the case of true 2n 
decay of ^^O than it was typically observed in the case of 
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FIG. 3. Convergence of a width F as a function of general- 
ized (hyperspherical) angular momentum K. Shown for the 
"moderate repulsion" case. Solid and dashed curves show cal- 
culations with Et ~ 150 and Et = 1.5 keV. For calculations 
above K — 20 the hyperspherical channels are taken into ac- 
count adiabatically. 



the 2p radioactivity (true 2p decay). For illustration how 
this mechanism is engaged we performed the "1/4 n-n 
FSI" calculation with the n-n potential multiplied by the 
0.25 factor. Though giving only 4 times larger width than 
the "no FSI" calculation at Et = 150 keV, this model at 
smaller decay energies rapidly sticks to the [s^]o system- 
atics providing dramatically larger widths. The triple- 
ridge correlation in the coordinate space pattern is now 
dissolved at much smaller radii [p ~ 80 fm, see Fig.^c)], 
compared to the "no FSI" case. The energy correlation 
here shows one strongly asymmetric peak indicating the 
[s^]/[p^] configuration mixing, which is consistent with 
the trend of the [s^] systematics. 

The full n-n FSI provides further drastic increase of 
the width (the dash-double-dotted curve in Fig. [2). The 
energy dependence is now completely following the [s^] 
trend. The coordinate correlation patterns indicate very 
fast (just above p ~ 30 fm) transition from the [<f] to 



TABLE I. Parameters of the Woods-Saxon potential in ■^*0-n 
channel for different model cases. Radial parameters are in 
fm and Vi are in MeV. 



Case 



Vs Vr, 



ro 



"1/4 n-n FSI", "Moderate repulsion" 70 70 3.5 0.75 

"no FSI" , "Strong repulsion" 120 120 5 1.2 

"No repulsion in p-wave" 70 3.5 0.75 
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FIG. 4. (Color online) Correlations in "^^O in different models: (a) "no FSI, no recoil", (b) "no FSI", (c) "1/4 n-n FSI" (d) 
"moderate repulsion". Each row at two left panels shows the coordinate space correlation density for ^^O WF vj((+) on two 
different scales. Two right panels show the energy and angular (momentum space) correlations in "T" and "Y" Jacobi systems. 
All cases are calculated with Et ~ 75 keV. 



the [s^] configuration, see Fig. lU^d). The energy corre- 
lation between neutrons is now only slightly asymmetric 
indicating a dominance of the [,s^] configuration at the 
asymptotics already at Ex = 75 keV. 



VIII. EFFECT OF THE OCCUPIED ORBITALS 

In our model we use the repulsion in s- and p-waves 
as a simple method to account for the effect of the occu- 
pied orbitals in ^^O in the absence of the experimental 
information about these interactions. To understand the 
importance of this aspect of the ^^0-n interaction for the 
^^O lifetime, we have performed calculations with very 
strong repulsive interactions in the s- and p-waves (the 
"strong repulsion" case) and without any nuclear inter- 
action in the p-wave (the "no repulsion in p-wave" case). 
The widths in these calculations are different by a factor 
of 50 (see Fig. [2|), while the lifetime energy systematics 



and correlations patterns are not varying on a significant 
level. We may consider the band between the latter two 
cases as the realistic prediction range defined by the un- 
certainty of this aspect of our model. It should also be 
noted that the attractive s-wave interaction may further 
increase the provided widths. However, at the moment 
there is no evidence for such a possibility in the available 
experimental data. 



IX. DISCUSSION 

Absence of Coulomb interaction in the case of long- 
living 2n emitters, compared to the 2p case, does not 
simplify, but makes the problem more complex in the 
sense of the sensitivity to several not well studied aspects 
of few-body dynamics. Nevertheless, we establish much 
narrower boundaries for the ^^O lifetime thus providing 
the conservative upper limit of ^^ 1 keV for the decay 



energy (assuming the lifetime from [13j). 

If the experimental ^^O lifetime is confirmed, the the- 
oretically derived decay energy is very small on nuclear 
scale: the energy of ^^O practically coincide with the 2n 
threshold. This would mean extreme "halo like" struc- 
ture of ^^O. Formally, the radius of the decaying system 
is infinite. Practically, for systems with radioactive life- 
times the radial characteristics are reliably saturated for 
integration in the subbarrier region. We have found val- 
ues around 5.5 fm for the "valence" neutron rms radius in 
^^O. Such values are typical for ^-'^Li which possesses the 
most extreme 2n halo known so far, which is connected 
mainly with the important [s^] component of the WF. In 
^^O the analogous radial properties are expected despite 
the [cP] component is dominating in the nuclear interior. 

The expected from experimental lifetime and our cal- 
culations decay energy of less than 1 keV underlines the 
importance of special experimental techniques aimed at 
studies of such phenomena. The direct measurement of 
the ^^O decay energy can provide an important cross 
check for possible discovery of the 2n radioactivity in [Ij] . 
However, such small decay energies are not accessible for 
existing experimental setups: what we know so far is just 
the fact that ^^O decays via 2n emission with intractably 
small decay energy. A method, based on the precision 
measurements of neutron angular correlations in reac- 
tions with relativistic secondary beams, could overcome 



this problem. This method was proposed in Ref. [14 | 
for future searches of 2n and 4n radioactivities and it is 
naturally suited for extreme low decay energies. Concep- 
tually analogous method for decays with proton emission 
was found to be very robust jl8i] . 



X. CONCLUSIONS 

The detailed studies of possible long-living (radioac- 
tive) true 2n emitters are performed for the first time by 
example of the ^^O system. The following main results 
should be emphasized. 

(i) The fine few-body effects play extremely important 
role in the decay dynamics of the true 2n emitters. The 
sensitivity of decay width to (a) configuration mixing due 
to core recoil, (b) subbarrier configuration mixing caused 
by n-n FSI, and (c) occupied orbitals effects far exceeds 
the corresponding effects in true 2p decays. Unexpect- 
edly, the lifetime systematics of ^^O with [cP] internal 
structure sticks to the typical [s^] behavior due to config- 
uration mixing, see Fig. [TJb). Studies of the correlation 
patterns provide deep insights in the 2n decay dynamics, 
(ii) The performed theoretical calculations provide much 
narrower limits for the lifetime vs. decay energy depen- 
dence of true 2n decay of s-d shell nuclei than those found 
in Ref. [ij|. These narrower limits also provide much 
more stringent limits on the decay energies at which the 
2n radioactivity may be found. 

(iii) We discuss the recently observed evidence for a 2n 
radioactivity of ^^0 (its lifetime is reported to be in a pi- 
cosecond range). Being confirmed such a lifetime means 
extremely low decay energy of ^^O (below ^ 1 keV) and 
calls for special experimental techniques for further stud- 
ies of this phenomenon. 
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